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Abstract 

Nonlinear dispersionless equations arise as the dispersionless limit of well know inte- 
grable hierarchies of equations or by construction, such as the system of hydrodynamic 
type. Some of these equations are integrable in the Hamiltonian sense and appear in the 
study of topological minimal models. In the first part of the review we will give a brief 
introduction to integrable models, mainly its Lax representation. Then, we will intro- 
duce the dispersionless limit and show some of our results concerning the two-component 
hyperbolic system of equations such as the polytropic gas and Born-Infeld equations. 



1. Introduction 



The study of integrable models or solvable nonlinear partial differential equations 
is an active area of research since the discovery of the inverse scattering method [1- 
3]. These models are in a sense universal since they show up in many areas of physics 
such as solid state, nonlinear optics, hydrodynamics, field theory just to name a few. 
Also, integrable models are linked to many areas of mathematics (see the chart in 



http : //www.ma.hw. ac .uk/ solitons/procs/bulloughl/bulloughl/bulloughl .html ) and have 



beautiful structures behind them. 

In this review we want to approach the dispersionless limit of some integrable models 
and describe some of our work on this subject [4-7]. This review is organized as follows: In 
Section 2 we review or at least introduce some basic facts on integrable models. We use the 
Korteweg-de Vries equation (KdV) as an example. In Section 3 we introduce the disper- 
sionless limit of an integrable model using the KdV equation to obtain the corresponding 
Riemann equation. Section 4 reviews our work with a special class of dispersionless sys- 
tems known as two-component hyperbolic systems. We show our results concerning the 
Hamiltonian structures for the Riemann equation [4] the dispersionless Lax representation 
for the poly tropic gas dynamics [5] and Born-Infeld equation [6]. Finally, in Section 5, we 
conclude with some problems that deserve further investigations. 

2. Integrable Models 
2.1 Solitons 

We are interested in nonlinear partial differential equations such as the sine-Gordon 
equation, nonlinear Schrodinger equation, Korteweg-de Vries equation (KdV), etc. These 
equations, as we will see, are very special since they are integrable. From now on we will 
illustrate the main results concerning integrability using the KdV equation. 

The KdV equation has as solution what is called today a soliton. We can trace the 
discovery of the soliton back to 1834 with the Scott Russell's experiment [8] to generate 
solitary waves in water, i.e., localized single entity waves. A modern version of his experi- 
ment is shown in Figure 1 (see jittp: //www.ma.hw. ac .uk/^chris/scott_russel .html| for an 



attempt to recreate Scott Russell's soliton). Scott Russell found that the volume V of 
water wave is equal to the volume of water displaced and that the speed c of the solitary 
wave is related with its amplitude a, depth of water h and acceleration of gravity g by 



c^=g{h + a) (2.1) 



This equation shows that higher waves travel faster. Attempts to obtain (2.1) theoretically 
were done by Boussinesq (1871) and Lord Rayleigh (1876) but an equation for u{x, t) in the 
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Figure 1. Generation of a solitary wave. 

small amplitude {h ^ a) and in the long wave regime {h ^ £) was deduced by Korteweg-de 
Vries in 1895 [9]. This is the now famous KdV equation 

~t~ ^XXX /„ „N 



where u{x, t) is the wave profile and ut = = ^ , . . . . 

The interest in the KdV equation (2.2) was resumed after studies of Fermi, Pasta 
and Ulam in 1955 [10] on numerical models of phonons in non-linear lattices, which are 
models closely related with the discretisation of the KdV equation. Motivated by these 
results, Zabusky and Kruskal in 1965 [11] studied numerically equations like (2.2) with 
periodic boundary conditions and were led to introduce the concept of "soliton" solutions. 
In 1967 Gardner, Greene, Kruskal and Miura [12] solved equation (2.2) exactly, introducing 
the "Inverse Scattering Transform Method" (ISTM), and were able to obtain its analytic 
expression. The so called 1-soliton and 2-soliton solutions of the KdV equation (2.2), for 
rapidly decreasing boundary conditions 

u{x,t) for a; — > ±oo , 



are 



u{x,t) =-c^sech^ (-c{x + c^t) ) 1— soliton 



/ N ,„3 + 4cosh(2a;-8t) +cosh(4x-64t) 

u(x,t) =12— — ^ — ^ — / ^ 2-soliton 

^ ^ {3 cosh X - 28t) + cos 3x - 36t }2 



(2.3) 



In Figure 2 we have pictures for the time evolution of the KdV soli- 
tons (2.3) (for some brief solitons movies see [http : //www . ma . hw . ac . uk/solitons| and 



http : //www . physics . otago . ac . nz/PhysicslOO/ simulat ions/Gamelan/ j ava/toda ) . The 1- 



soliton solution in Figure 2 is the solitary wave obtained in the Scott Russell's experiment. 
Observe that as the time evolves the wave keeps its form. For the 2-soliton solution in 



KdV 1-soliton 



KdV 2-soliton 



Figure 2. Time evolution for the solitons of the KdV equation. 

Figure 2, since the taller the soliton the faster it moves, the two solitons will interact nonlin- 
ear ly when they meet. But, the amazing fact is that the two solitons will almost keep their 
initial form after interaction, there will be only a shift in their positions. This particle-like 
character and ability to retain its identity after interactions is what characterize a soliton 
solution of a nonlinear equation such as the KdV one. 



2.2 Inverse Scattering 



The next breakthrough in the soliton thread came in 1968 with the Lax [13] discov- 
ery about the meaning of the ISTM. His observation is that the KdV equation has the 
representation 



dL 
'dt 



[S,L] 



where 



L =d^ + -u 



B=4d^ + -{du + ud) 
2 



(2.4) 



(2.5) 



are operators. Here d = ^ satisfies df = + fd. We call L the Lax operator and in 
some sense we can find a Lax representation such as (2.4) for any integrable system. In 
this way, starting from (2.4), we can apply the ISTM for other nonlinear equations. 
We can write the following eigenvalue problem for the Lax operator L 



Lip = -\ip 



(2.6) 



It is easy to see that since L evolves in time as (2.4) we have = 0, i.e., the eigenvalue 
problem is isospectral. For the KdV equation (2.6) assumes the form 



^+(-Ma:,t)+A)V = 



(2.7) 



which is the time-independent Schrodinger equation and where t is a parameter (not the 
time in the Schrodinger equation). Now we can obtain a solution u{x,t) as follows: For 
some given initial condition u{x^ 0) we solve (2.7) and obtain the scattering data S{t = 0), 
since u satisfies the KdV equation we can obtain the scattering data for any t, so from 
S{t) we use the inverse scattering (as we usually do in quantum mechanics) to find the 
"potential" u{x,t) from the scattering data S{t). This is the ISTM routine and the main 
steps are illustrated in the diagram bellow. 
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Inverse Scattering Transform Method 



2.3 Hamiltonian Systems 

In 1970 Gardner [14] showed that the KdV equation is a Hamiltonian integrable 
system. Then, Faddeev and Zakharov in 1971 [15] were able to interpret the ISTM as a 
change of variables to the action angle variables. In fact, the representation of integrable 
models as integrable Hamiltonian systems is the starting point to the "Quantum Inverse 
Scattering Method" . Before we see how the KdV equation can be expressed in Hamiltonian 
form let us review the symplectic formalism for Hamiltonian systems. A Hamiltonian 
system is described by a phase space qi,Pi, with i = 1, . . . , N, and a Hamiltonian function 



H(pi,qi). The equations of motion are then given by the Hamilton's equations 

dH 



Pi 



dpi 
dH 

dqi 



(2.8) 



Alternatively, we can describe a Hamiltonian system using Poisson brackets, for the dy- 
namical variables A{q,p) and B{q,p), defined by 



{A,B} 



OA dB OA dB 



(2.9) 



dqi dpi dpi dqi 

which is skew-symmetric and satisfies the Jacobi identity. The variables of phase space 
satisfy the canonical relations {qi,qj} = {pi,Pj} = and {qi,Pj} = Sij. The Hamilton's 
equations (2.8) assume the form 

qi ={qi,H} 
Pi ={Pi,H} 

Putting the variables q^ and Pi in an 2N dimension column z the equations (2.8) assume 
the form 

/ qi\ ( 9/ dqi \ 



(2.10) 



d_ 
di 



qN 
pi 

\pN/ 
= Z 



/ 
-/ 

"=7 



d/dqN 
d/dpi 

\d/dpN) 



H 



(2.11) 



or 



and even in a more compact form as 



z = JVH 



(2.12) 



(2.13) 



This is the symplectic formalism for Hamiltonian systems. The Poisson brackets can be 
written as 

(2.14) 



{A,B}= {VA^^ J {VB^ 



where J"'^ = —J^"' and ^ (y^J'^^d^J^^ + cyclic) = 0. The canonical relations are given by 
{z,z} = J and (2.10) by 



(2.15) 



We can perform some generalizations, aUowing J to depend on 2;, J{z), and going 
from a discret sympletic space, of dimension 2A'", to the continuum where we have now a 
field u{x,t) instead of z{t). Then, we have the following "dictionary" 



z{t) 
H{z) 

VH 

J{z) skew-symmetric matrix 

z = JVH 

{z,z} = J{z) 



{A,B} = (VA^ j(vB 



u{x, t) 

H[u] functional 

SH ^ . , , . . 

functional derivative 

ou 

T>{u) skew-adjoint operator 
u = V- 



Su 



{u{x),u{x')} — VS{x — x') 



If there is a J ^ we say that we are in a symplectic manifold, otherwise we are in a more 



general situation of a Poisson manifold. Note that the functional derivative 



5H\ 



Su 



as 



5H[u{x)] _ ^.^ H[u{x) + e6{x - y)] - H[u{x)] 
Su{y) e^o e 



is defined 
(2.16) 



which for H[u] = u{x) yields 

SH[u(x)] 
5u{y) 

and for H[u] — J dx h{x, u, Ux, Uxx, ■ ■ •) 

5H[u(x)] fd d d 



5u{y) 



5{x - y) 



du dx du^ dx^ dui 



h 



where the right hand side is just the Euler-Lagrange operator acting on h. 

Now, let us return to the KdV equation (2.2) and observe that it can be rewritten as 



'U't — UUx ~\~ Uxxx 

d (\ 2 

-u -|- Un 



dx\2 



(2.17) 



Introducing the Hamiltonian 



H2= I dx( ^u' - \ul ) (2.18) 



we see that = ^tt^ + u^x and = 0. The operator 



= I- (2.19) 

is skew-adjoint and satisfies the Jacobi identity. So, (2.17) can be written in Hamiltonian 
form as 

Ut = V^—^ = {u{x),H2}i (2.20) 
ou 

where 

{u{x),u{y)]i=Vi5{x-y) (2.21) 

and we are omitting the exphcit dependence on t. 

Besides (2.18) the KdV equation (2.2) has an infinite number of conserved charges 

Hq = [ dxu 
Hi = I dx -v? 



/dx -\ 
2 



H2= I dx i - \ui 



3! 2 



=jdx 



= I dx { -u^ — Zuux + ^w^j. 



(2.22) 



, , , 1 5 ^ 2 2 36 2 2 

H4^= I dx [ -u - 6m + —uu^^ - -^u^xx 



and it can be shown that these charges are in involution, i.e., 

{if„,if^}i = (2.23) 

making the KdV equation integrable in Lioville's sense. 

In 1978 Magri [16] discovered that equations like KdV have a second Hamiltonian 
structure. The operator 

1 ( d d 



is skew-adjoint and satisfies Jacobi identity, and the KdV equation can be written in the 
alternative Hamiltonian form 

{u{x),Hi}2 (2.25) 

V2S{x - y) (2.26) 

These charges (2.22) are also in involution with respect to this second Hamiltonian struc- 
ture 

{Hn,Hm}2 = ^ (2.27) 

We say that the KdV equation is a bi-Hamiltonian system. In general we say that a system 
is bi-hamiltonian if there are Hamiltonian operators X'l and 'D2 which are compatible, i.e., 
such that Pi, 1^2 and \iDi + A2l^2 satisfy the Jacobi identity. It can be shown [16] that 
if a system is bi-Hamiltonian it is integrable in Lioville's sense. 

Starting with the works of Gel'fand and Dickey in 1975 [17], Adler in 1979 [18] and 
many others, algebraic developments started to take place. The key role played by the 
Lax operator L, in obtaining the conserved charges Hn, the Hamiltonian structures, the 
hierarchy of equations that share was then revealed. In the next sections we will 
introduce and apply some of these techniques in the dispersionless situation. 

3. Dispersionless Limit 

We have seen that solitons preserve their shape and speed after collision. The soliton 
solution has a nondispersive nature. This is so not because dispersion effects are absent 
but because there is a compensation by the nonlinearities of the system. Let us look at 
the KdV equation (2.2) more closelly. If we eliminate the nonlinear term in (2.2) we get 
the linear dispersive equation 

Ut = Uxxx (3.1) 

which admits the solution 

u{x, t) dk A(/c)e^('=^-^('=)*) (3.2) 

This is a pure dispersive solution. In Figure 3 we see that a initial configuration at t = 
will disperse as time goes on. Eliminating the dispersive term we get the pure nonlinear 
equation 

Ut = uux (3.3) 

It can be easily checked by substitution that 



ou 

where 

{u{x),u{y)}2 = 



u{x, t) = f{x — ut) 



(3.4) 



with / arbitrary, satisfies (3.3). From this solution we conclude that the velocity of a 
point of the wave, with constant amplitude u, is proportional to its amplitude leading to 
the "breaking" of the wave, as shown in Figure 3. The wave also develops discontinuities 
(indicated by the vertical dashed line in Figure 3) in its evolution. The "miracle" of the 
soliton solution is due to a balance between the dispersion and the breaking of the wave, 
both phenonema placed together lead to the wave profile to propagate without changing 
its shape. 




Figure 3. The balance effects of dispersion and breaking in a soliton. 

Equation (3.3) is called the dispersionless KdV or Riemann equation [19]. The inter- 
esting fact is that this equation is a integrable Hamiltonian system. We will return to study 
this equation in the next section but for the moment let us analyse how we get dispersionless 
equations. Dispersionless equations can be obtained by construction or as a quasi-classical 
limit of integrable ones [20] . In the latter case we make the scaling ^ a^, |^ ~^ '^§^ 
and take the limit a — > 0. For the KdV equation (2.2) (we will change the constant factors 
on it, so instead of (2.5) we have L = + u and B = + j{du + ud)) 

4ut = Uxxx + ^uUj: =^ AaUf = a^Vfixx + QcxUxU 

JJ. Ut = \uUx 

KdV ^ 



Riemann 



This is like the WKB aproximation in quantum mechanics and we will use it as our guideline 
[20]. 

Dispersionless integrable systems were introduced by Lebedev and Manin [21] and 
Zakharov [22], and although interesting on their own started to appear recently in de- 
velopments in low-dimensional quantum field theory. It has been shown that there is a 
connection between 2-dimensional field theories and integrable equations of hydrodynam- 
ical type [23-25] (which are dispersionless systems). In 2-dimensional topological field 
theories [26] we are interested in calculating, from the partition function 



Zm^ J [#] e-^[^l (3.5) 



the correlation functions 

{4'a{x)(t)i3{y) • • •)m = (0a0/3 " " ■) M (3.6) 

which depend only on the topology of the manifold M. The 2-point and 3-point correlation 
functions are given respectively by [27] 

{4'a4'f3) = Vap = nondegenerate constant 

d^F{t) (3-7) 

{(pa(Pl3(pf) = Cap^it) = Q^aQ^HQi-y a, (3, 'y ^ 1, 2, . . . , U 

where t = (t-*^, t^, . . . , t") are the coupling constants and F{t) is the free energy. The 
correlations (3.7) define a commutative and associative algebra (with an identity) 

Ca o = clfje^ (3.8) 

with Cq; defining a basis for the algebra. The associativity of the algebra, (cq, o e^j) o = 
Ca o (e/3 o e^), gives 

d^F{t) d'Fjt) _ d^Fjt) d'Fjt) 

I Qj.'^ Qj-A Qj-y/ Qj.'^ Qj./^ Qj-A ' Qj-rv- Qj-A Qj./;. \^'^ ) 



These are the Witten-Dijkgraaf-Vcrlinde-Verlinde (WDVV) equations [26,27] and can be 
identified with equations of hydrodynamic type. So, solutions of hydrodynamic equation 
can be identified with particular solutions of the topological field theory [25] . 

4. Two-Component Hyperbolic Systems 



In a series of papers [28-31] Nutku and collaborators started to study dispersionless 
systems of equations that are in-between the simple Riemann equation [19] and the more 
general equations of hydrodynamic type [24] 



Riemann 



Two Coinponciit Hyperbolic 



Equations of Hydrodynamic Type 



ui = V^^^ui i,j = l, 



,n 



In Figure 4 we can find a chart with the main equations, of the two-component hyperbohc 
system type, studied on these papers. 

A wealth of results concerning the integrability of these systems were revealled. In- 
finitely many conservation laws and multi-Hamiltonian structures were obtained. In this 
section we will be interested in reproduce some of these results from an algebraic point of 
view. In order to achive this goal we must understand the Lax representation for these 
systems. 

4.1 Riemann Equation 

The Riemann equation 



(4.1) 

is the prototype for the hyperbolic systems. We address the following question: Is there a 
Lax representation for (4.1)? Yes, and we can obtain it performing the semiclassical limit 
[20] explained in Section 3. So, if the KdV equation goes to the Riemann equation (4.1) 
in the semiclassical limit, the Lax operator L — + u and B — + j{du + ud) goes to 
the polynomials in the variable p 



E = p'^ + u 

3 3 
M = p'^ + -up 
2 



(4.2) 



and the Lax representation (2.4) goes to 



(4.3) 



called dispcrsionless Lax representation (note the resemblance when we pass from quantum 
to classical mechanics doing d ^ p and [ , ]—*■{,})• Here 



{A{x,p),B{x,p)} 



dA dB dB dA 



dp dx dp dx 



(4.4) 



is the dispcrsionless Poisson bracket [21,22]. So, if we substitute (4.2) in (4.3) we get (4.1). 
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H{u,v) = ^ + l 
i 
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Figure 4. Two-component hyperbolic equations. 



From now on we will apply some of the techniques described in [17] and [18] in a very- 
informal way, since we want to give only a flavor of how the "machinery" works. 

Let us calculate the square root of E in (4.2). So, we write the Laurent polynomial 



£;V2 =p_^(^^_^aip ^ +a2P ^ + a3P ^ + 



(4.5) 



and from E = E^/'^E^/'^ we obtain ao, ai, 02, as, . . ., or equivalently, we perform a series 
expansion for p — > cxd 



u 



1/2 



Now we calculate E^/'^ = E^/^E, E^/^ = E^/^E^ and so on 



^3/2 _^ 3 ^ _^ ^u'^P~^ + ■■■ 
2 8 



5 ±u o 

z o ib 



15 



,3^-1 



(4.7) 



The set of general Laurent polynomial A = J^t^oo ^iP^ gives rise to an associative algebra 
g = {A}. This algebra can be written as a direct sum g = (B Q-, where g^ — {^+} ^md 
g- = {A-} with A+ = Yli>o '^iP^ ^- = Si<o '^iP^^ respectively. We can recognize M 
in (4.2) as 

M = {E^^^)+ (4.8) 
In fact, from (4.3) we are motivated to write 



^ = {{E^/^U,E} 
dE 



dt 

dE 



Ut = Ua 



3 

15 



(4.9) 



and we have a hierarchy of equations. We call it dispersionless KdV (or Riemann) hierarchy 
and we write 

dE 2fe+i 

— = {{E—)+,E}, A; = 0,l,2,3,... (4.10) 



treating tt as a function of A; + 1 variables 

u = u{x,to,ti,t2,...) (4.11) 
For each t}. we have what is called a flow and it can be shown that they commute 

(4.12) 



dtidtk dtkdti 

consequently, the whole set of equations (4.9) is integrable since, as we have already pointed 
out, the Riemann equation is an integrable Hamiltonian system (all the equations in (4.9) 
share the same set of conserved charges). 

The Riemann equation can be put in the form 

3 3 

ut = -uux = -^{u'^)x (4.13) 

It follows that the quantity H oc J dxu^ is conserved. In fact / dx vT' are conserved as we 
can show explicitly. These conserved charges can also be obtained from E. Let be A any 
general Laurent polynomial 

A=--- + a_ip-^ + ■■■ (4.14) 
following [18] we introduce the Adler's trace as 

Tr^ = J dxResA = j dxa-i (4.15) 

which satisfies the usual relation Ti:AB = TrBA. From (4.6) and (4.7) we see that 



16 J 



and we have 



(4.16) 



whit Hn ~ 0. From a Hamiltonian point of view the Riemann equation is a quadri- 
Hamiltonian system [30]. There are Hamiltonian operators T>i, T>2, T>^ which are com- 
patible and another Hamiltonian operator £ which is compatible only with T>i. We can 
write 



Ut = Vi 



SHr, 

5u 



5u 



3^ '^^i 
4 ^ Su 



—8 — - 
8 5u 



(4.17) 



where 

Hi ^ J dxu, Vi = 2d 

H3 = - dxu^ , V2 = ud -\- du 



H^ = - I dxu\ Vs = u^d + du'' ^^-^^"^ 



Hq = — dxu^, £ = d-^d^d 

Hamiltonian structures can also be obtained from the Lax operator L {E in the disper- 
sionless case) . They are the symplectic structures of Kostant-KiriUov [32] on the orbits of 
the coadjoint representation of Lie groups [18,33]. For dispersionless equations the corre- 
sponding Lie algebra is given by the associative algebra of Laurent polynomials endowed 
with the bracket (4.4). For the KdV equation the Lie algebra is given by the algebra of 
the pseudo-differential operators with the usual commutator. Following this scheme the 
Hamiltonian structures "Di, 'D2, can be derived (see [4] for details) while we were not 
able to obtain £ from this scheme. 

4.2 Polytropic Gas Equation 

We will try to apply the results of the last Section to some others dispersionless 
equations, such as the ones in the chart of Figure 4. The polytropic gas dynamics equation 



Ut + UUx + v'^ ''vx =0 
vt + {uv)x =0 



7 > 2 (4.19) 



was studied from a Hamiltonian point of view in [30]. In (4.19) u is the velocity of the 
fluid, V is its density, / = v'^~'' and is related to the pressure (/(f) = ^-^) and 7 is the 
ratio of specific heats (we call an ideal gas polytropic if the specific heats are constant over 
a large range of temperature). 

The first step will be to derive a Lax representation for (4.19). We get a hint if we 
consider 7 = 2 in (4.19). In this case we have the shallow water equation [19] also known as 
the irrotational Benney equation [34] . Even though we do not know the dispersive system 
which originates (4.19) for any 7 we do know it for the case 7 = 2. This is the dispersive 
shallow water [35] equation, also called the two boson equation in field theory 

§^=(2J, + Jo2 + J^)' 

/ (4.20) 
^=(2JoJi + JO' 



This equation has the following nonstandard Lax representation [36,37] 



L =a - Jo + d-^Ji 



dL 
'dt 



(4.21) 



where (L^)>i stands for the purely nonnegative (without terms) part of the polynomial 
in p and Jq oc tt, Ji oc t; are the two bosons fields. Now, if we perform the semiclassical 
limit and do the appropriate identifications (4.20) yields (4.19) for 7 = 2 and from (4.21) 
we get the following dispersionless Lax representation 



L =p + u + vp 



(4.22) 



For any 7 we can use (4.22) as an ansatz to obtain the dispersionless Lax representation 
for (4.19) and it reads [5] 




(4.23) 



In [30] two sets of conserved charges were derived for (4.19) when 7 7^ 2. So, if (4.23) is 
really the correct Lax pair it must somehow provide both sets accordingly to the algebraic 
scheme described in the last section. In fact, since L has singularities in p — and p — 00 
we can expand in powers of p in the two following ways 



Lt-1 = p^l + 



7-1 
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So, the first set of charges follows from 




n = 0,1,2,3, 



where the first densities are 
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and 
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which are the first set of charges obtained in [30]. The second set follows from 



(4.27) 




n = 0,1,2,3, 
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and the first densities are 

Ho =(7- l)-^v 
Hi =(7 — 1) T-i —uv 
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where 
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is the second set of charges obtained in [30] 
1 

In (4.23) Lt-1 was expanded in p = oo, a expansion around p = provides a second 
consistent dispersionless Lax equation 



dL 
'dt 



(4.31) 



which yields (with the proper rescahng) the equations 



■7-3 



Vt = Ua 



(4.32) 



From the chart in Figure 4 we recognize this equations as the polytropic elastic media 
equation. 



4.3 Born-Infeld Equation 

With the Lax representation for the polytropic gas, obtained in the last section, we 
can get a Lax representation for the Born-Infeld equation given in the chart of Figure 4 




(4.33) 



In (4.33) the Born-Infeld equation is expressed in the so called null coordinates version 
[31]. If we perform the transformation 



w = 0a 



V = 



(4.34) 



VI + 4>x4>t 

we obtain the Born-Infeld equation written as a second-order equation in null coordinates 



4>'i4>tt + 4>'i4>xx - (4 + 24>x4>t)4>xt = o 



(4.35) 



A Lax representation for (4.33) can be obtained as follows [6]. In the first place if we do 
the change of variables 

u = -(u^ +v'^) 

^ 1 (4.36) 



V = -uv 
2 



called Verosky transformation [31], we will end up with the equation 



Ut + UUx + — = 
Vt + iuv)a; = 



(4.37) 



known as the Chaplygin gas. In view of this it would be desirable to first obtain a Lax 
description of the Chaplygin gas like equations 



=0 , a> 1 



(4.38) 



Vt + {uv)a; =0 

This is indeed possible if we set 7 — > —a, a > 1 in (4.23), so (4.38) can be obtained from 

(4.39) 
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where is expanded around p = and (L^Tt)^^ is the polynomial in p that produces 

consistent equations, instead of the purely nonnegative polynomial used in (4.23). For 
a = 1 the Lax operator 



L=p-^- (4 + 4) +^P^ 



dL 
'dt 



reproduces (4.33). Again, conserved charges follows from 

= Tr L"^^^ = J dxHn n = 0,l,2,3. 



and the first Born-Infeld charges are 
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and these are exactly the charges derived in [31]. Another set is obtained from 

-- j dxHn n = 0,l,2,3,... (4.43) 




and the first ones are 
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(4.44) 



This is a new set of conserved charges, for the Born-Infeld equation (4.33), not found 
previously in [31]. 

5. Conclusions 

We believe, from the results of the Section 4, that the study of dispersionless systems 
via a Lax representation is worthwhile. So, the search for a dispersionless Lax representa- 
tion for the equations in the upper part of the chart in Figure 4 is being pursued. Also, 
the derivation of the multi-Hamiltonian structures of these systems, as described in [28- 
31], is under investigation following the coadjoint orbit method [32,33]. Another question 
that comes to mind is the dispersive generalization of these equations. Attempts in this 
direction can be fovmd in [38]. 

Some topological equations are also related with the systems discussed here. For 
instance, the hyperbolic Monge- Ampere equation 

UttUa^a. - {Ut,f = -1 (5.1) 

may be related with the Born-Infeld equation as follows. If we perform the change of 
variables 

a — Ux 

(5.2) 

b = Ut 

the Monge- Ampere equation can be written as a first order system 

at = h 

hi - 1 (5.3) 



h 



a. 



and this equation can be related to the Chaplygin gas equation (4.37) through the following 
change of variables 

u = 

ax (5.4) 

V = ttx 

Thus, we can give a Lax description for the hyperbolic Monge- Ampere equation through the 
Lax representation derived in Section 4.3. Finally, the Witten-Dijkgraaf-Verlinde-Verlinde 
(WDVV) equations (3.9), for n = 3, with 

F{t\t\ t') = ^{t')H' + h\t^f + f{t\ t') (5.5) 

where = x and = t, yields the third order Monge- Ampere equation 

fttt = fxxt ~ fxxxfxtt (5-6) 

This equation is a bi-Hamiltonian system and has a matrix Lax representation. It is then 
possible to generate a whole set of nonlocal charges much like the nonlinear sigma model 
(details are given in [7]). It is likely that a dispersionless sort of Lax representation for 
(5.6) may exist. 
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